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Abstract 

A new class of entanglement witnesses (EWs) called reduction type entanglement 
witnesses is introduced, which can detect some multi-qudit entangeled states including 
PPT ones with Hilbert space of dimension (g) d 2 (g) ... (g) d n . The novelty of this work 
comes from the fact that the feasible regions turn out to be convex polygons, hence the 
manipulation of these EWs reduces to linear programming which can be solved exactly 
by using simplex method. The decomposability and non-decomposability of these EWs 
are studied and it is shown that it has a close connection with eigenvalues and optimality 
of EWs. Also using the Jamiolkowski isomorphism, the corresponding possible positive 
maps, including the generalized reduction maps of Ref. [21], are obtained. 
Keywords: Entanglement Witness, Multi-qudit, Optimality, Linear Program- 
ming, Feasible Region . 

PACs Index: 03.65.Ud 



Multi-Qudit EW 



3 



1 Introduction 

Quantum entangled states lie at the heart of the rapidly developing field of quantum informa- 
tion science, which encompasses important potential applications such as quantum communi- 
cation, quantum computation (TJ [2j However, the fundamental nature of entangled states 
has tantalized physicists since the earliest days of quantum mechanics, and even today is by no 
means fully understood. One of the most basic problems is that how can one tell if a quantum 
state is entangled?, and how entangled is it still after some noisy quantum process. 
Here we will deal with the states of finite dimensional Hilbert space 7i = Ti.^ <8> ... <E> T~td n - A 
density matrix p is called unentangled or separable if it can be written as a convex combination 
of pure product states as I7) = \a^) ® ... ® \a^). If no such convex linear combination exists 
for a given p, then the state is called entangled. Although, in the case of pure states of bipartite 
systems it is to check whether a given state is, or is not entangled, the question is yet an open 
problem in the case of mixed states (see the recent good reviews jU El El E])- There have been 
many efforts in recent years to analyze the separability and quantitative character of quantum 
entanglement. The Bell inequalities satisfied by a separable system give the first necessary 
condition for separability [H]. In 1996, Peres made an important step towards proving that, 
for a separable state, the partial transposition with respect to one subsystem of a bipartite 
density matrix is positive, p TA > [9J. By establishing a close connection between positive 
map theory and separability, Horodecki et al. promptly showed that this is sufficient condi- 
tion for separability for bipartite systems of 2 (g> 2 and 2 eg) 3 |10| . Regarding the quantitative 
character of entanglement, Wootters succeeded in giving an elegant formula to compute the 
entanglement of formation [H] of 2® 2 mixtures, thus giving also a separability criterion [T2] . 
An alternative method to detect entanglement is to construct so-called entanglement witnesses 
[TOl [TBI [T4l [T5] and positive maps [16] . Entanglement witnesses are physical observables that 
can detect the presence of entanglement. Recently many attempts have been made to use the 
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convex optimization methods as a robust tools in most quantum information areas including 
construction of EWs PH DSl EE HH H5] . 

The main motivation for the present investigation is the remarkable fact that there is no evi- 
dence that when and where the manipulation of the EWs reduces to LP which can be solved 
exactly. For example, the first author and his collaborators have found certain examples of 
Bell diagonal EWs and generalized d x d Choi entanglement witnesses [HI [15], by approxi- 
mate LP method. Indeed in most cases determining feasible regions needs to use numerical 
calculation and consequently the problem is solved approximately. Therefore to figure out 
the problems which can be easily solved exactly in a optimal way and also can be generalized 
to an arbitrary number of particles with different Hilbert spaces is the main purpose of this 
paper. To this aim we introduce a new class of (d 1 ® d 2 (g) ... Cg) d n )-multi-qudit EWs called 
reduction type entanglement witnesses(REWs) and show that the computational difficulty in 
solving such problems reduces to LP which can be solved by simplex method[20j. In these cases 
the feasible regions are simplexes and the minimum value of optimization problem is achieved 
on their apexes. On the other hand we show that the EWs corresponding to hyperplanes 
surrounding feasible regions are optimal. Another advantage of this work is the fact that all 
REWs can be written in terms of some positive operators and optimal EWs. Also in most 
cases the decomposability of REWs is rather determined, where two of their eigenvalues plays 
an important role in this issue. Another important consequences of such EWs are the positive 
maps including the generalized reduction map [2T] . which can be obtained from these EWs 
or their tensor product, via Jamiolkowski isomorphism. Finally a class of entangled density 
matrix including PPTES are provided which can be detected by such REWs and consequently 
they determine the non-decomposability of REWs. At the end we carry out in detail some 
particular cases such as multi-qubit and 2 <g> 3 <S> 4 and generalized reduction EW. 
The paper organized as follows: In section 2 we give a brief review of EWs. In section 3 we 
explain the general scheme of linear programming. In section 4, we introduce a new class of 
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EWs which can be put in the realm of LP, since their feasible regions are convex polygons 
(indeed simplexes) which can be exactly determined. In section 5, we show that all EWs cor- 
responding to hyperplanes surrounding the feasible regions are optimal. Section 6 is devoted 
to some interesting examples such as: multi-qubit REWs, (d<g>d...®d) multi-qudit REWs and 
2® 3® 4 REW. Section 7 deals with two examples of LP type, where the first one can be solved 
exactly by the prescription of this paper. In section 8 we introduce some entangled and PPT 
states which can be detected by REWs and the decomposability of REWs is discussed. In 
section 9 by using Jamiolkowski isomorphism, the relation between REWs and positive maps 
is explained. The paper is ended with a brief conclusion and three appendices. 



2 Entanglement witness 

As mentioned above in the introduction, one of the pragmatic approach to detect entanglement 
is to construct entanglement witnesses. Let us briefly recall what these operators are. 

Definition 1. A Hermitian operator W is called an entanglement witness detecting the en- 
tangled state p e if Tr(yVp e ) < and Tr(Wp s ) > for all separable state p s . 

So, if we have a state p and we measure Tr(Wp) < 0, we can be sure that p is entangled. 
This definition has a clear geometrical meaning. The expectation value of an observable 
depends linearly on the state. Thus, the set of states where Tr(Wp) = holds is a hyperplane 
in the set of all states, cutting this set into two parts. In the part with Tr(Wp) > lies the 
set of all separable states, the other part ( with TriWp) < 0) is the set of state detected by 
W. From this geometrical interpretation it follows that all entangled states can be detected 
by witness. Indeed for each entangled state p e there exist an entanglement witness detecting 

it pm. 

Definition 2. An EW is decomposable (d-EW) iff there exists operators V, Qi with 

W = V + Ql A + Ql B + ... + Q T N Z V,Qi>0 (2.1) 
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where partial transpose taken with respect to some subsystems and it is non- decomposable if it 
can not be put in this form [19]. 

Only non-decomposable EWs can detect PPT entangled states that is those density ma- 
trices which have positive partial transposition with respect to each subsystems [22J. 



3 Manipulating EWs by LP method 

This section deals with basic definitions of linear programming(LP) and general scheme to 
construct EWs by an exact LP method. 

To this aim first we consider a Hermitian operator W with some negative eigenvalues 

W = Y,a i a i (3.2) 

i 

where o~i are positive operators, with < Tr(aip s ) < 1 for every separable states p s . Note 
that, the condition < Tr(aip s ) < 1 is not always required. It is used here only to simplify 
analyzing the problem and pave the way to generalize the prescription to multi-qudits with 
arbitrary higher dimensions as will be discussed in the following. 

The minimum value of T = Tr(Wp s ) achieves for pure product state, since every mixed p s can 
be written as a convex combination of pure product states (due to the convexity of separable 
region) as p s = YliiPi\li){li\ with Pi > and J2iPi = 1' hence we have 

Tr(p s W) = J2Pi Tr 0^h)(^\) > C ™in (3-3) 

i 

with C min = minjF := Tr(yV\j)(-?\) \j) G D prod , 

where D pro d denotes the set of product states. Thus we need to find the pure product state 
\lmin) which minimize Tr(W|7)(7|). Now, as the pure product state I7) varies, the map 
defined by Pi = Tr(<jj|7)(7|) maps the set D prod into a region inside the hypercube defined by 
< Pi < 1 i = 1,2, ...,N. This is not all circumstances. An important and difficult task 
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is to find the convex region (called feasible region) inside this hypercube which comes from 
JF W := TriWps) as p s varies on D sep , where D sep denotes the set of separable states. Here 
in this work we are interested in the EWs with the feasible regions of simplexes (or at most 
convex polygons)types, such that the manipulating these EWs amounts to 

minimize T w := Tr(Wp s ) = '^^a i P i (3.4) 

% 

N 

subject to ^^(QjPj — di) > j = 1, 2, ... 

i=i 

where Qj ,di i,j = 1,2,. ..are parameters of hyper-planes surrounding the feasible regions. 
Therefore, the corresponding boundary points of feasible region will minimize exactly JF W , 
thus the required task reduces to LP problem which can be solved by the well-known simplex 
method [20]. 



4 (e?i ® g?2 ® ... ® ^ ra )multi-qudit reduction type EWs 

In this section we consider n Particles with arbitrary dimensions. Without loss of generality 
particles can be arranged so that d\ < d 2 < ... < d n . The discussion of some special cases 
is postponed to Section EJ We introduce and parameterize the multi-qudit reduction-type 
entangled witnesses (REWs) labeled by subscript R as 

w i n) = E h s °s + di & 2 ,...jvv..o><vv..oi + E «> ( 4 - 5 ) 

S^N' S^N' 

where N' = {2, ...,n} , with b % = b x , a' = a' and a s , a' s defined as 
di-l 

a s = J2 N) (*l ® °f ] ® - ® °t ] with °t ] = 

and 

< = EE - E M&l ® of ® ... ® q£»> (4.7) 

U=0 12=0 i n =0 



if fceS 
/ if jfe £ 5 



(4.6) 
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with 



O 



'(*) 



\h)(h\s ilih 

KfcX^IIlLfi^iC 1 -* 



o 



l»fcXifc|(l-*i*ii) 
respectively and 



if fceS 

,,,,, w if fc ^ S 1 , |5|<n-3 and i 2 , 
if k S, \S\ = n — 2 and z 2 , 
otherwise 

di— 1 



(4i 



i=0 



is maximally entangled state. In this notation we have a' = a[. Obviously for multi-qubit 
system none of a' exists. The number of P' s depends on the dimensions of particles d^s and 
it can take one of the following values 



m 



E 



|S|=n-d C |s , | 1 if di — d 2 



En— 2 (~m—l 
|ol_n M 



. = d n = d 

,151=0^151 ~ (2 m_1 — 1) if di = ... = d m < d m+ i < ... < d n 
We introduce the new more convenient parameters a 1 = b 1 = b% , a s = b s + E s , Cs b s ,, instead 
of 6 s 's. In order to turn the observable (14.51) to an EW, we need to choose its parameters in 
such a way that it becomes a non-positive operator with positive expectation values in any 
pure product state (7) = \a^') (g> ... ® |cr n '). 

Now it is the time to reduce the problem to LP one. Io order to determine the feasible region, 
we need to know the apexes, namely the extremum points, to construct the hyperplanes 
surrounding the feasible regions. Suppose that (7) be a pure product state with |ckW) = 
(a( fc ),a[ fc ),...,aW_ i ) T and let 



di-l 



P s :=Tr(a s \ 1 )( 1 \) = J2\^ 1) ^ ) -^ 



(n)|2 



a 



W if ke{l}uS 



i=0 



1 if k£{i}us 



P 2,...,n ■= d x |(Vv..ol7>F 



di-l 

E 

i=0 



P^:=Tr(a'| 7 )(7l), 
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which all lie in the interval [0, 1] (see appendix I). The number of Ps is 2 n_1 — 1 (cardinality 
of the power set of N' excepted {0}). 
The extremum points or apexes consist of 

a. origin: P s = 0, P' s = VS* g N' which corresponds to the following choice of pure 
product state 

= (100 ... 0) T 
\a^) = (010 ... Of k G N' 

b. p s = i, p s > = i vs" g s, p N , = o, p s » = o vs" gN'\s, p' s = o vs g iV' 

which can be reached by choosing the following pure product states 

= (100 ... Of k G {1}US, 

= (0 1 ... Of k(£{l}US. 

Obviously if Ps = 1, then for all S" g S, Ps> = 1, thus for 5 = iV', we get the following 
important apex 

c. P s = 1, P^ = VS g TV 7 , 
which can be obtained by choosing 

\ a ( k )) = (100 ... 0) T \/k G {l}ur 



d. Pg = 1, P s > — 1 VS" g S g TV , the others are zero. 

The last category arise from pure product state with the components of the form = a, 
1, if k G S and aj^i = 1, if /c ^ 5* , i.e., 

= (100 ... Of fc G {1}US, 

= (0 ... ^ ... Of k {1} U 5. 

fc-l'th 
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Regarding the above consideration, we are now ready to state the feasible region. 
To this aim we first prove that, 2 n_1 + m extremum points obtained above form the apexes of 
(pn-\ _|_ m )_ s i m plex in Euclidean space of dimension N = 2 n_1 + m — 1. To this purpose we 
consider the convex hull of theses points, i.e. draw N + 1 hyperplanes passing through each 
combination of iV points out of N + 1 ones( C^ +1 = N + 1). Now, we get a bounded region 
formed from their intersection which is the required feasible region of W^ ra - ) and it is obviously 
a (2 n_1 + m)- simplex. It is strait forward to show that the feasible region can be obtained by 
taking the expectation values of 



(4.9) 

in pure product states (see appendix II). Now, in order that Wj^ to be an EW, the following 
expectation values 

f(P 2 ,P 3 ,...,P[,P>,...):=Tr(W^\ 1 )( 1 \) = £ h, Ps ■ b., ,r-i ,. • Y. a 's P 's ( 4 - 10 ) 

Sgv' S<^N> 

must be positive. So our task is to solve the following LP problem 



minimize 



E b S Ps + K.,A-,n+ J] < P' s 



Ps-Ps + H (-1) |5|+|5,| P 5 , > 
subject to { sc * s ' ~ (4.11) 

V P' s > 

Putting the coordinates of apexes in Eq. (14.101) yields all a s > and a s + a' s > 0. As we stated 
in previous section, all Ps and P$ lie in the closed interval [0, 1]. Now, T is a linear function 
of Ps and P' s and if we require it to be positive on the apexes ( which are extremum points), 
then it will be positive in the whole feasible region. 

At the end we need to know all eigenvalues of Wj^ which consist of a s ,a s + a' s ,cu 1 = a N , — 
b N , and u 2 = dia N , — (d\ — 1)oj 1 . Since a s ,a s + a' s > 0, then one of the remaining eigenvalues: 
uj 1 and u 2 must be negative to guarantee W^ n ^ to be an EW. 
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5 Optimality of {S) W^ 

After determining the feasible regions, one needs to know whether the boundary of EWs is 
formed by optimal EWS. An EW is optimal if a positive operator P is subtracted from that 
then it will be no longer an EW [T3] . 

Note that the EWs corresponding to hyperplanes surrounding feasible regions of Wj™ are 
optimal since they cover the simplex feasible region in an optimal way (see appendix III). 
Thus, the structure of the optimal EWs characterizes the boundary of REWs VVj^. 

In fact, from the results of this Section it will become clear that we can restrict ourselves to 
the structure of the optimal EWs corresponding to hyperplanes surrounding feasible regions. 
In other words optimal EWs ^VV^ are tangent to the boundary between separable and non- 
separable states. 

Another advantage of ^VV Q y is that one can rewrite the in terms of positive operators 
a s' lVv..o)(VV..ol anc ^ some optimal EWs, i.e., 

=E a S (5) <" ) + rf 1 ^l^o...o)(^,.,l + EK + <)< S^N' (5.12) 

S S 

Therefore for positive u 2 , the REWs can be decomposed as 

)V(«) = ^(s)qVw +7) (5 . 13 ) 
s 

where ( S 'Q := ^VV^/ (JV ' XS) and V is positive operator and in this case Wj^ can not detect 
PPT entangled states( non-separable density matrices with positive partial transpose with 
respect to all particles). 



6 Some special cases of 

R 

In this section we discuss some special cases of REWs to enlighten the subject. 
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6.1 Multi-qubit reduction type EWs 

It is important both theoretically and experimentally to study multi-qubit entanglement and 
to provide EWs to verify that in a given multi-qubit state, entanglement is really present. The 
Eq. (14. 5 p for system of n-qubits reduces to 



(6.14) 



As mentioned before, the dimension of qubit-systems does not allow presence of a' s . 

The number of P's is 2 n_1 — 1 , whereas the number of apexes is 2 n ~ 1 . Putting the 
coordinates of apexes in Eq. fl4.10l) again indicates that all a s > 0. The feasible region is 2 n ~ 1 - 
simplex of dimension 2 n_1 — 1 surrounded by hypersurfaces defined by (14. lip with all P"s 
eliminated. 



6.1.1 Three-qubit (n=3) 



The first nontrivial example of REWs for multi-qubit system is three qubit REW 
= a i J 8 + ( a 2 - a i) a 2 + ( a 3 - a iK + 2 Ks +a 1 -a 2 - a 3 )\ip 000 )(ip 000 \ 



(6.15) 



with eigenvalues > 0, u> 1 = a 2 + a z — a 1 and u 2 = 2 a 23 — u x . In order to obtain the 

feasible region, we need the coordinates of apexes which are 



It) 



It) 



It) 



It) 



/A /,\ / \ 











(P 2 = 0,P 3 = 0,P 2 , 3 = 0) 



(P 2 = l,P 3 = 0,P 2i3 = 0) 



(P 2 = 0,P 3 = 1,P 2 , 3 = 0) 



(p 2 = i,p 3 = i,p 2i3 = i; 
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subject to 



As we see the number of Ps is three, while the number of optimal points is four. Therefore, 
there are four hyperplanes surrounding the feasible region, where the hyperplanes pass through 
each combination of three points out of four ones, (see Fig.l). Thus the problem can be reduced 
to 

minimize Tr(W^\ r y)('y\) 

1 + P 2 ,3 " P2 ~ P3 > 0, 
P2 ~ ^2,3 > 0, 

(6.16) 

P3 ~ P2,3 > 0, 
P2,3 > 0, 

where the above given inequalities follow rather easily by taking the expectation value of the 
following optimal EWs 

3 Wgl=a 3 ( v 3 -2\rl> m )U> fm \) (6-17) 
!>vg) = Qi (/ 8 _ a2 _ az + 2 I ^000)^000 1) 

in pure product states, and as usual the optimal EWs can be written as partial transpose of 
the following positive operators 

1 wgr 23 =a 1 (|100) + |011))((100| + (011|), 

2 <f 3 = 0,(1001) + |110»((001| + (110|), (6.18) 
3 wg T 2 = a 3 (|010) + |101»«010| + (101|), 

respectively. Now the EW Wj 3 ) can be written in terms of % W opt as 



WW = a^Wopt + a 2 2 W opt + a 3 3 W opt + 2 co 2 ^ 000 )^ 
= aWQ T(23) + a 2 ( 2 )Q T(3) + a^Q T{2) + 2 uj 2 V 



(6.19) 



As mentioned above in section 4, for cu 2 > the EW W^ 3 -* becomes decomposable and can not 
detect PPT entangled states. 
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6.2 n-qudit (d® d® ... ® d) 

For n particles with the same dimensions the extra terms a' will appear in EWs, provided 
that the requirement \S\ < n — 3 is met. Then the Eg. (14.51) becomes 

K l) =Y. b s°s + d k jvv..oX</v.J + E a >' s ( 6 - 2 °) 

S^N> S<=JV',|S|<n-3 

where 



< = Et!,^.Ni)^l«o? ®...®£3£", of 



We discuss below the most simple case: an REW consisting of just two qudits with the same 
dimension d, that is 

Where for a 2 = it reduces to the well known reduction EW ( the term "reduction-type EWs" 
for general Wj^ is inspired from this particular case). In this case we have only P<i which 
can takes values between and 1. So the feasible region is just the line segment < P-i < 1. 
Putting and 1 in T{P>i) = a 1 + (a 2 — ajP 2 yields a 1 > and a 2 > respectively. The 
eigenvalues are a 1 and 2 a 2 — a 1 where the second one must be negative to ensure detecting 
some entangled states. 

6.3 Three particles (2 <8> 3 <8> 4) 

As a particular example of REWs with different dimension let us discuss three particles with 
2 ® 3 <8> 4 dimensions 

= a 1 / 8 + (a 2 -a 1 )a 2 + (a 3 -a 1 )a 3 +2 (a 2 ^+a 1 -a 2 -a 3 )|^ 000 )(^ ool+«+ a 2 CT 2+ a 3 (T 3 ( 6 - 21 ) 
where we have 



a[ = |012)(012| + |013)(013| + |021)(021| + |022)(022| + |023)(023| 
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|102)(102| + 1 103) <103| + |120)(120| + |122)(122| + |123)(123| 

a' 2 = |002)(002| + |003)(003| + |112)(112| + |113)(113| 

<j' 3 = |020)(020| + |121)(121| 

Here all possible a' s (S = {2}, {3} and 0) can be appear. In this case feasible region lies in a 
space of dimension six. The coordinates of apexes and relevant pure product states (7) are 



(P2, P3, P23, P[i P'21 P3) 



|7> 



|7> 



|7> 



V / 



V / 



|7) = 



V / 



\ 



1 





v / 



1 

v / 



(A 



v / 



<8> 



1 



/ \ 
1 





v°y 






(0,0,0,0,0,0) 



(1,0,0,0,0,0) 



(0,1,0,0,0,0) 



(1,1,1,0,0,0) 
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It) 



It) 



It) 



v / 



V / 




1 

v / 





V / 
i 

V / 





1 

V / 
/o\ 



1 

V / 






(0,0,0,1,0,0) 



(1,0,0,0,1,0) 



(0,1,0,0,0,1) 



Again choosing all combinations of six apexes out of seven ones, one can find the boundary of 
feasible region as 

P{,P^P^P 2>3 >0 

P 2 - P 2i3 - P' 2 > 
P3 - P2,3 ~ Pi > 
1 + P2,3 -P2-P3-P{>0 

where the EWs corresponding to these hyperplanes are 



2 W opt = a 2 (a 2 -2 | ip 000 ){ip 000 1 - cr' 2 ), 
3 W opt = a 3 ( a 3 - 2 | ^000)^000 1 - °3)> 
^opt = a^I^ - a 2 - a 3 + 2 |^ 000 )(^ 000 | - oi). 



(6.22) 
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Taking the partial transposition of t W op t, i — 1,2,3 with respect to {2, 3} \ {i} yields 

l Wo$ = ai(lioo) + |011))((100| + (0111), 

2 WS = 0,(1001) + |110»«001| + (110|), (6.23) 
3 WS = a 3 (|010) + |101))((010| + (101|), 

respectively. Evidently these EWs are optimal, since these have been written as the partial 
transposition of pure maximally entangled states. 



7 Bell-diagonal EWs by LP methods 

Recently multi-qubit Bell decomposable entangled witnesses (BDEWs) [14] have been intro- 
duced as 

W = a . . . 1-0- ■ ■ )(V>. ■ - I (7.24) 

i 1 i2---in— 0,1 

where |'0j li2 in ) (^i — 2, % — 1, 2..., n) are n-qubit maximally entangled orthonormal states, i.e., 

IVW..J = ® ^) i2 ® - ® K) in l^oo...o) (7-25) 

where cr^ and a z are usual Pauli matrices. 

In general it is hard to find the BDEWs with feasible region of simplex type or even polygon 
type, namely those which can manipulated by LP method. Here we give two examples which 
are both set in LP problem, where only one of them (the first example) can be solved exactly 
by the prescription of this paper. 

The first example is EW of the form 

Wi = al 2 u + 2{b - a)|V 00 ...o>(Vv..ol + 2 ( c ~ a )l^oo...oi)(^oo...oil + (d-a)a (7.26) 

where 

a = I 2 n - {(|VV... 110 )(VW..l 1 ol + IVV..lo)(VV..lol) + (l^011...l)(^011...ll + IV> 1 l...ll)(^ll...ul) 
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00. ..01 00. ..01 I + IV>10...0l)(V>10...0ll)} 

The eigenvalues of Wi are a, 26 — a, 2c — a, d. This BDEW is similar to the one introduced 
in [TJ], where the extra term a is added to optimize the EWs corresponding to the boundary 
plane of feasible region. Suppose that P 00 = 2 \(ip 00 |7)| 2 , -P 00 01 = 2 \(ip 00 01 |7)| 2 and 
P = Tr(cr|7)(7|). Then the pure product states which produce the apexes are 









'°1 


|7> = 








|7> = 

< 








|7> = 


>:< 


® 




|7) = 

















«„.,, =0, />„,.„= 0,P = 0) 



v / 








v / 



(p =0 P =0 p = v 

\ 00. ..0 w > 00. ..01 w > / 



(p =1 p 

\ 00. ..0 ' 00. ..01 



fp = P 

V 00. ..0 w ' 0( 



0,P = 0) 



1,P = 0) 



and consequently these yield the following hyperplanes surrounding the feasible region (see 
Fig.2) 



P P P > 

00. ..0 ' 00. ..01 J — w 



1 - P -P -P>0 

00. ..0 00. ..01 — w 



The positivity of the last constraint come from the positivity of the expectation value of the 
following optimal EW 



W opt = /-2|Vv.. )(Vv..ol - 2 IVw..oi)(Vw..oil -° 



in pure product states I7), since it can be written as the partial transpose of a positive operator 
with respect to the fist particle as 



wS = 2 (IV'ii...io)(V'i 1 ..,ol + l^i...ii)(^i...iil)- 
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Now the remaining task is to solve the following LP problem 

minimize a + (b — a)P 00 + (c — a)-P 00 ... 01 + (d — a)P 
1 - P - P - P > 

00. ..0 1 00. ..01 1 — w , . 

subject to (7 .27 ) 

P P P > 

00.. .0 > x 00. ..01 ' — w 

Thus, above problem is reduced to LP and can be solved by simplex method. Putting the 
apexes in Eq. ( 17.261) we deduce that a, b, c, d should be positive. Now, the operator Wi fulfills 
the properties of EWs if at least one of its eigenvalues is negative, namely 2b — a < or 
2c- a < 0. 

The second example which sets in LP problem is 

W 2 = al 2 n + 2 n -\b - aMo...o>(Vv..ol + - a)|W..o>(Vv..ol 

which can not be solved by the prescription of this paper. Its feasible region can be determined 
by Lagrangian multiplier method, as it is discussed in [14]. Let P00...0 = 2 | (Vv..o It) | 2 >-Plo...o = 
2 |(^ 10 |7)| 2 . There the problem reduces to 

minimize a + 2 n ~ 2 (b - a)P 00 + 2 n_2 (c - a)P 10 



subject to 



2"- 



£z - P,o.., + (1 - ^)^oo...o > (7-28) 
P P > 

± 00. ..0 ) 10. ..0 — w 



These constraints can not be reach by partial transposition approach, and the feasible region 
is estimated by convex hall of apexes (see Fig. 3). 

8 Detecting some entangled states by Wj^ 

This section is devoted to some entangled states which can be detected by general and 
three-qubit REW . First we consider some Bell states. All of the Bell states ) can 
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be detected by since we have 

therefore for u 2 < one can detect \ip 00 ) and for cu 1 < the others modulated Bell states 
can be detected by Wj™- 1 . On the other hand, imposing some constraints on operator 

„ = BTr{p ) + Ddl | B ^ ~ E 1^ o...o)(^ 0...0I) + ^ dilV 4 o...o>(^ 00...0I I , (8-29) 

one can get a density matrix which can be detected by W^, where p s denotes the separable 
state inside the parenthesis on the righthand side. The positivity of p i0 constrains B and 
D to be positive and in order to detect both p Q0 and p. ,z 7^ 0, we must have 

B Tr(W^p s ) + Dd 1 u 2 <0, (8.30) 

B Tr(W^p s ) + Dd 1 u 1 <0, (8.31) 

respectively. Because of the positivity of D, Eq. (18.301) is satisfied if u 2 *C and Eq. (l8.3ip 
is satisfied if co>i < 0. Now, we can proceed our discussion further to detect PPT entangled 
states which is useful for determining non-decomposable region. Here we do not discuss the 
decomposability and non-decomposability issues in detail, since it needs the other opportunity 
and comes elsewhere. p i is PPT states with respect to any subsystems of the particles, if 
B — D be positive, so the Eq. (18.301) yields 

B —d\ uo 

1< — < Vt— => uj 2 <zu 8.32 

~D Tr{wP Pa ) 



where 



w:=~Tr(W^p.) (8.33) 



The above requirement makes p a class of PPT entangled state which can be detected by 
Wj™* 1 . But to detect the other p i0 we must have ui\ < w which is impossible. This is in 
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agreement with the discussion made in section [5j Now, combining thus obtained results with 
those of Eq. (15.121) which implies that Wj™ is decomposable provided that u 2 is positive, one 
can rather determine the decomposability and non-decomposability of REWs (see Fig.4). 
Furthermore one can construct some entangled states which can be detected by particular 
REWs. As an example consider entangled density matrices 

P = 4j g j 2jp (£ ^ + 2 WoooX^ooo I) 

' {Ba,+2D\^ Qm )^ J) 



AB + 2D 

which can be detected by three-qubit REW (16.151) . with some constraints. The positivity of 
these states implies that B, B + 2D > and the positivity of p TA , p' T2 is achieved if B ± D > 0. 
In order that, Tr(W^ p) to be negative we should have B(a 2 3 + a 2 ) + Dlu 2 < 0, where we have 
two possibilities: for D > we have uj 2 < and for D < we have uj 2 > 0. 



9 Positive maps 



As it is shown in [23] , there is a close connection between the positive maps and the entangle- 
ment witnesses, i.e., the Jamiolkowski isomorphism 



di < di 



i — 1, n 



£{p) = Tr ( i i3) ..., 2 n-i) I wSJ"::^,^^'®-^^"-^) 



where 



d\— ld.2 — 1 (in — 1 



y/did 2 ...d n 



EE-Ei i( i 1, 4 2, )i4 3, 4 1, )...i€"- 1, C" ) ) 

11=0 12=0 i n =0 



(9.34) 



(9.35) 



(9.36) 



is the maximally entangled state in TC^ <8> W^/ ® ■■■®'H^ <S>7"^™' ) . Hence the Jamiolkowski 
isomorphism is a one-to-one mapping between the set of trace preserving quantum operations 



e : <> » < » ... » «2T» — ««? ««' ... ® «>> 



(9.37) 
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and di x d! X ... x d n x d' n EWs if d{ < d[ for i = 1, 2, n. 

Now, using the Jamiolkowski isomorphism (19.351) . we try to construct the positive maps con- 
nected with REWs. Evidently tensor product of some EWs is also an EW in higher dimension. 



;(l,2,...,2n) K ^ _ pw _ + ._ _ ^(1) ^ ^(2) ^ ^ ^(2n-l) ^ ^(2n) 

depending on possible partition 



To be more precise let W££^£j, be an EW acting on H ( d \> <g> 7^7 <g> ... <g> ft^ n j <g> ft^f j then 



n = n 1 + n 2 + ••• + w m , > 1 (9.38) 
one can construct an EW by tensor product of REWs as 

yy(l,2,...,2n) = yy(l,2,...,2ni) ^y(2ni+l,...,2ni+2n 2 ) „ „ yy(2n-2n m +l,...,2n) 

cfi,^,...,^,^ di,d' 1 ,...,d ni ,d^ 1 d„ 1+ i,d' ni+1 ,...,d ni+n2 ,d' ni+n2 d n _ nm + i,^_ nm+1 ,...,d n ,d^ ' 

(9.39) 

then using Jamiolkowski isomorphism (19.351) . one can obtain the corresponding positive map. 
For instance, considering the tensor product of n REWs(corresponding to the partition n = 
1 + 1 + .. . + 1) 

n 

w (l,2,...,2n) _ X-X w (2fc-l,2fc) 
yV d 1 ,d' 1 ,-,dn,d' n - \&) yV d k ,d' k 
k=l 

with 

yy^" 1 ' 2 ^ _ a (2fc-l,2fc) j + ^ a (2fc-l,2fc) _ a (2fe-l,2fe)^|^(2fe-l,2fe)^^(2fc-l,2fc)| + a /(2fe-l,2fe) ./(2fe-l,2fe) 



acting on Hilbert space ® 7i^, with a'^ 2k 1,2k ^ given in Eq. (14. 7ft . we get the following 

positive map 



where 



S^{l,3,...,2n-1} 



r s = II (^-a^), 

JGV\S 
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and (ji...j\s\) is the ordered of S. Choosing all a 2 = a[ = and a l = 1 yields the generalized 
reduction map introduced in [21]. As an example for n = 2 one can easily verify that 

+a^\a^ - a^)Tr s (p) ® J« + (af- 2 ) - a^)(a^ - af^)p 

d' 2 -l d'. 2 -l 

flM^TrWjf ® E l*><*l + (°?' 2) " «S 1,2) )< (3 ' 4) rr 3 (/>) ® E |*><*| 

k=d,2 k=d.2 

d[-l d[-l 
+a ;(^) (a (3,4) _ a (3, 4)) J- l^^l g, Tri(p) + a ' i (l,2) a (3 I 4) Tr(p) J- l^^l f (4) 

fc=cii k=di 

d{-l 4-1 
+G /(l,2) a /(3,4) Tr(p ) ^ l^^l g, ^ |^| 

k=di k=d2 

Taking ap) = af- 4 ) = af > 2 ) = af > 4 ) = and af- 2 ) = aj 3 - 4 ) = 1 yields 

E®(p) = Tr(p)l£> ® /£> - /£> ® Tr^p) - Tr 3 (p) ® /£> + p (9.40) 
and choosing a* 4,2 ) = a^ 3,4 -* = and a'- 1,2 -' = a^ 3,4 ) = — a'S 1 ' 2 ' = — a'( 3 ' 4 ) = 1 we reach the new 

° 2 2 11 1 1 

reduction positive map for different dimensions 



d' 2 -l d' 2 -l 



£( 2 \p)=Tr(p)lV®l£ ) -I^®Tr 1 (p)-Tr^ \k)(k\+Tr 3 (p)®J2 \ k )( k \ 

k=d2 k=d2 



d[-l d[-l d{-l 4-1 

+ \ k )( k \ ® t ^i(p) - E ® 4? + Tr (p) E i*x*i ® E \ k ^ k \ 

k=d\ k=d\ k=di k=d2 

These examples show that one can construct more positive maps by making tensor product of 
EWs in arbitrary way provided that the dimensionality condition di < d\ is satisfied. 



10 Conclusion 



The generalized reduction type entanglement witnesses with simplex feasible regions are in- 
troduced, where the EWs corresponding to hyperplanes surrounding the feasible regions are 
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optimal. These REWs are of types that their manipulation is reduced to LP problem and can 
be solved exactly by using the simplex method. As it shown above, the REWs are decompos- 
able in cases if their second eigenvalue, namely uj 2 becomes positive while for negative values of 
oo 2) the decomposability or non-decomposability of REWs is still open for debate. Also various 
other interesting issues remain unsolved, such as keeping the REWs in realm of LP problems 
despite of adding some other operators or entangled states to them. 

Appendix I: 

Proof of the inequalities: < Ps,P's,P2,...,n < 1- 

In this appendix we prove that all Ps , P' s and p2,..., n take the values between and 1. The 
inequalities < Ps, P' s < 1 can be easily concluded from the following ones 



0<Tr(a s | 7 )( 7 |)<nEl^ fe) | 2 = 1 ' 



fc=l i=0 



0<TrK| 7 )( 7 |)<nEK (fe) | 2 = 1 - 



k=l i=0 

For i-2,...,n) the Cauchy-Schwartz inequality implies that 



di-l L 

P 2 ,...,n := di KVw..ol7>| 2 = E a ! 1)a ! 2) - a ! n) =\(a (1) \p)\ 2 < |||a (1) )iril|/?)|| 2 , 



where 





finally using the following inequality 



di-l 



n dfr — 1 




i=0 



k=2 i=0 
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one can conclude that < p2,..., n < 1- 

Appendix II: 

Proof of the feasible region of ( 14.111) . 

In order to obtain the feasible region of (14.111) . we need to evaluate the expectation value of 
optimal EWs 



\ S^S'^N' / 

(1-1) 

in pure product states |t) (t I where =^ VV^. Now by taking the partial transpose of 

( s )yV (r ^ with respect to (N' \ S) we have 



opt 

(ii-l 



^W( p f™=a s ^|^f)(^f| (1.2) 

i+3 



where 

AS) 



1*8°) := l«?> ® l«!?> ® - ® + ® |/?<?>) ® ... ® |/?W>, (1-3) 

I I*) if k E IU S ... I | j> if kelUS 

I | j) if fc^lUS Mi) if k <£ 1 U 5 

Noting that all of these operators are positive definite and using the relation 



Tr{ (5) W ( P " } (JV ' XS) |7)(7l ) =Tr{ (5) W ( J(|7> (l\) T ^ ) > 
yields all feasible regions which are simplexes. 



Appendix III: 

Proof of the optimality of ^W^. 

Here in this appendix we try to prove that witness is optimal, to this aim we give 

the proof for the special case ( 2 )yy( n ) since the proof of general case is rather similar to this 
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particular one. According to the Reference [33], the EW ^VV^, is optimal if for all positive 
operator P and e > 0, the following new Hermitian operator 

W new = (1 + e)< 2 >W<2 - sP (1-5) 

is not anymore an EW. Suppose that there is a positive operator P and e > such that 
W new = ^Wj™ — eP is yet an EW. Let the positive operator P be the pure projection 
operator since an arbitrary positive operator can be written as sum of pure states 

with positive coefficients as P = J2i -M^X^il- 

Now, one should note that the expectation value of the operator ^'W™ in pure product 
states 1 7) will vanish if they satisfy the following equation 

AiB* + AjB* = (1-6) 

with 



Ai = (a 1 ) i (ot a ) i , Bj = (a 3 ) j (a 4 ) j ...(a n ) j . 

But, it is straightforward to see that, for A iy Bj G 1R, the pure state will be similar to 

one of the with i ^ j with 

\Mf) ■■= \i) ® |»> ® Ij) ® - ® Ij) + Ij) ® Ij) ® |*> ® - ® |*> i,j = o,...,di - 1, (1-7) 

concluding that an arbitrary P has the form: P = a ij 

l*y> Wl with ^ > 0. 

Finally, substituting Eq. fll— €3[) in the following expression 
Tr(P\-y)(j\) = ^a^-IA^- + A,-^ 2 = J2 a ij 

ij ij 

and choosing £>j's such that -BJI^ 7^ -Bj-B* yields a i:; - = and consequently one can conclude 
that P = 0. 



B„ 



BjB*\ 
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Figure Captions 

Figure-1: 3-simplex displaying the feasible region of three-qubit REW. 
Figure-2: 3-simplex displaying the feasible region of multi-qubit BDEW Wi. 
Figure-3: Convex polygon displaying the boundaries of the feasible region for multi-qubit 
BDEW W 2 . 

Figure-4: Decomposable and non-decomposable regions of REWs: for uo 2 > the REWs 
are decomposable, for uo 2 < w the REWs are non-decomposable and in dashed region, the 
decomposability or non-decomposability of REWs is still open for debate. 
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